Numerical analysis of pressure oscillations over axisymmetric spiked blunt bodies at Mach 6.80

The pressure oscillations over a forward facing spike attached to an axisymmetric blunt body are simulated by solving time-dependent compressible Navier–Stokes equations. The governing fluid flow equations are discretized in spatial coordinates employing a finite volume approach which reduces the equations to semidiscretized ordinary differential equations. Temporal integration is performed using the two-stage Runge–Kutta time stepping scheme. A global time step is used to obtain a time-accurate numerical solution. The numerical computation is carried out for a freestream Mach number of 6.80 and for spike length to hemispherical diameter ratios of 0.5, 1.0 and 2.0. The flow features around the spiked blunt body are characterized by a conical shock wave emanating from the spike tip, a region of separated flow in front of the hemispherical cap, and the resulting reattachment shock wave. Comparisons of the numerical results are made with the available experimental results, such as schlieren pictures and the surface pressure distribution along the spiked blunt body. They are found to be in good agreement. Spectral analysis of the computed pressure oscillations are performed employing fast Fourier transforms. The surface pressure oscillations over the spike and phase plots exhibit a behaviour analogous to that of the Van der Pol equation for a self-sustained oscillatory flow. Received 28 February 2001 / Accepted 17 January 2002     

A Field Monte Carlo Formulation for Calculating the Probability Density Function of a Single Scalar in a Turbulent Flow

The probability density function (PDF) formulation of one scalar field undergoing diffusion, turbulent convection and chemical reaction is restated in terms of stochastic fields. These fields are smooth in space as they have a length scale similar to that of the PDF. Their evolution is described by a set of stochastic partial differential equations, which are solved using a finite volume scheme with a stochastic source term. The application of this methodology to a particular flow is shown first for a linear source term, with exact analytical solution for the mean and standard deviation, and then for a nonlinear reaction.     

沟槽面与光滑面湍流边界层特性比较

应用ＬＤＶ技术对沟槽面及光滑面湍流边界层流速和湍流度分布进行了精细的测量，实验结果表明：与光滑面相比，沟槽面湍流边界层时均流速分布对数公式中具有较大的积分常数Ｃ值，且沟槽面湍流度最大值较小，但其出现的位置距壁面较远。另外，偏斜因子及平坦因子的分布特性区别不大，但是在近壁区内沟槽面结果波动较大。     

Reynolds stress, mean velocity, and dynamic static pressure measurement by a four-hole pressure probe

An improved version of the four-hole directional pressure probe, or Cobra probe, is described, in which the frequency response has been extended to 1.5 kHz. The probe measures all three orthogonal mean and turbulent velocity components at a point in the flow field. The probe also resolves the local mean and turbulent components of static pressure, allowing moments between the fluctuating velocity components and pressure to be determined. The techniques developed to allow the improved frequency response and the use of the probe in turbulent, developed pipe flow (a calibration flow) are described. Also given are the turbulent pressure-velocity correlations, which show a high degree of anticorrelation for one velocity component.     

湍流边界层中马蹄形涡的形成方式

本文使用一种新的流动显示方法——激光片光运动法和几种实验技巧对湍流边界层中的马蹄形涡进行了观测,发现并描述了其形成的四种方式:二次不稳定式、组合式、变形式和突发式。对这四种马蹄形涡的形成及发展进行了研究和比较。实验结果表明,这些马蹄形涡在尺度、运动速度和变形上是有差别的。     

防护林带：湍流的数学模型与计算机模拟

虽然防护林用于减小风速、控制热量和水汽传递及污染物扩散、 改善气候与环境、增加作物产量等已经有几百年了, 但直到近几十年, 人们才开始系统地研究防护林空气动力学的遮蔽机制.在本综述中,我 们考察了绕防护林带的流动与湍流控制机制,最新的模型与数值模拟 研究情况;通过数值模拟与实验数据的比较,来了解防护林带结构与防 风效果之间的关系;讨论数值分析如何及为什么能够得到所需要的结果. 本文将从多孔隙防护林带流动基本方程组的推导开始,讨论数值模型及 模拟过程,对附体与分离流动进行预测;分析了遮蔽机制与动量交换;对 风向、防护林密度、宽度和三维性对流动与湍流的影响作了系统的论 述.还对热流和土壤水分蒸发的新模型及数值模拟作了简述.最后,我们 对网络工作站、群和高性能分布式并行计算机及其对防护林带模型预 报能力的提高作了讨论.     

Numerical modeling of the disturbances of the separated flow in a rounded compression corner

Numerical modeling of the time-dependent supersonic flow over a compression corner with different roundness radii is performed on the basis of the solution of the two-dimensional Navier-Stokes equations in the regimes corresponding to local boundary layer separation. The development of unstable disturbances generated by local periodic injection/suction in the preseparated boundary layer is calculated. The results are compared with those of similar calculations for a flat plate. It is shown that the natural oscillations of the boundary-layer second mode stabilize in the separation zone and grow intensely downstream of the reattachment point. The acoustic modes excited within a separation bubble are studied using numerical calculations and an asymptotic analysis.     

Numerical Experiments in Scalar Diffusion under Stochastic Convection

We study scalar diffusion, both from Eulerian and Lagrangian perspectives, advected by two dimensional flows. Emphasis is devoted to the problem of scalar diffusion under a synthetic turbulent flow. We present numerical and analytical results for the turbulent diffusion coefficient either from the influence of the turbulent synthetic field and a periodic array of eddies. Preliminary results concerning Lagrangian dispersion are also reported.     

Vorticity shedding from a lentil-shaped body at large incidence in uniform flow

The 2D flow about a lentil-shaped body at high incidence in uniform stream is numerically simulated for high Reynolds numbers. The separation points are fixed on the body edges and the vorticity released through a Kutta-like condition is spatially discretized by point vortices. Their Lagrangianinviscid dynamics is described by using both aconformal mapping and aboundary elements method. The present results are in good agreement with previous calculations but still not satisfactory when compared with the experiments.The aim of the present work is to investigate the effect of different Kutta conditions and the effect of different amalgamation procedures in the far wake on global quantities such as theStrouhal number or the total shed circulation. Other remarkable effects induced either by the action of viscosity or by the rising of secondary 3D motion, may not be analyzed by the present simplified model.

---

Sommario Si studia il flusso attorno ad una lente ad elevata incidenza nelle ipotesi di elevato numero di Reynolds e flusso bidimensionale. I punti di separazione sono fissati sulle estremità della lente e la produzione di vorticità da questi è realizzata con opportune condizioni di Kutta. La dinamica dei vortici è descritta con un modello Lagrangiano non viscoso. I risultati, ottenuti sia con l'uso di una trasformazione conforme che con un metodo agli elementi di contorno, sono in ottimo accordo con i risultati di precedenti simulazioni numeriche, ma non sono ancora soddisfacenti se confrontati con gli esperimenti.Si analizza in termini di grandezze globali del flusso, quali il numero di Strouhal e la circolazione totale rilasciata, l'effetto di due diverse condizioni di Kutta e dell'uso di procedure di compattazione della scia lontana. Si conclude che l'uso di diverse condizioni di Kutta produce risultati pressoché identici mentre l'utilizzo di procedure di compattazione della scia può alterare in modo significativo il numero di Strouhal del flusso. Rimangono da analizzare altre possibili cause del disaccordo con le sperimentazioni quali la diffusione viscosa e la probabile presenza di moti tridimensionali.

     

Cloud dispersion models

A new technique has been developed to compute mean and fluctuating concentrations in complex turbulent flows. An initial distribution of material is discretised into any small clouds which are advected by a combination of the mean flow and large scale turbulence. The turbulence can be simulated either by Kinematic Simulation or by a stochastic model for the motion of each cloud centroid. The clouds also diffuse relative to their centroids; the statistics for this are obtained from a separate calculation of the growth of individual clouds in small scale turbulence, generated by Kinematic Simulation. The ensemble of discrete clouds is periodically rediscretised, to limit the size of the small clouds and prevent overlapping. The model is illustrated with simulations of dispersion in uniform flow and in a coastal flow, and the results are compared with analytic, steady state solutions where available.